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A typographical error was discovered in the computer code used to generate the results presented in the paper. The 
corrected results are displayed in Table I, which also incorporates new results, obtained using the recently improved 
post-Newtonian calculations of Tanaka, Tagoshi, and Sasaki 

The quantity computed in Table I is 

■f^-. ^ _ ^/A^lactual ^ mK)\ 

G^ 

the ratio of the actual signal-to-noise ratio obtained when matched filtering a gravitational-wave signal h with a 
O \ post-Newtonian template /i„ (n denotes the order of the approximation), to the maximum possible value. The precise 
O ■ definition of the inner product (-j-), and all other technical details, are given in the paper. 

Ph ' The reference signal h is constructed from exact expressions for both dE/dt, the rate at which energy is lost to 
gravitational waves, and dE/df, which provides the relation between orbital energy and gravitational- wave frequency. 
dE/dt is obtained by numerically integrating the Teukolsky equation for a small-mass object in circular motion around 
a Schwarzschild black hole. dE/df is obtained by integrating the geodesic equations for circular orbits. 

On the other hand, the template /i„ is constructed from post-Newtonian approximations to both dE/dt and dE/df, 
^ [ accurate through order w", where v = {ttM f)^^^ (with M the black- hole mass) is the orbital velocity. The expressions 
■ used are llj 

^ ! = ^(^) = 1 - 3-7113i;2 -f 12.566w^ - 4.9285w^ - 38.293^;^ -I- (115.73 - 16.305 \tlv)v^ - lOl.bW 
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+ (-117.50 + 52.743 In t))u^ + (719.13 - 204.89 lnu)v'' + (-1216.9 + 116.64 Inw)^^ 
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where {dE/dt)^, {dE/df)M are the leading-order, Newtonian expressions. The first form for Q{v) is exact. The 
second form (fully expanded) is a Taylor expansion about v = Q. The third form (partially expanded) forces Q{v) 
to vanish dX v — 6"^^^, the innermost stable circular orbit. The paper's original calculations used post-Newtonian 
expansions accurate through 0{v^). 

In calculating 7^„, the mass and arrival time parameters are assumed to be the same in h and /i„. In the original 
calculations, the initial-phase parameters were also assumed to be matched. Here we proceed differently: the overlap 
integral {h\hn) is maximized over all possible values of the difference in initial-phase parameters. While this obviously 
leads to larger results for 7?.„, the improvement is not significant: at most 2% when 71 > 6. 

The paper's conclusions are drastically affected by the coding error. 

First, the values obtained for TZn are now much larger than what was obtained before, and the Table shows that TZn 
nicely converges to unity as the order of the post-Newtonian approximation increases. In particular, the Table shows 
that a post-Newtonian template accurate through order would reproduce more than 95% of the signal-to-noise 
ratio. (System B is an exception.) 

Second, there is no longer any evidence for the observation that an accurate knowledge of the location of the 
innermost stable circular orbit would significantly improve the performance of the post-Newtonian templates. Indeed, 
in most cases the numbers displayed in the fourth column of Table I, calculated using the partially expanded form for 
Q{v), are not significantly larger than those displayed in the second column. (And in some cases they are significantly 
lower.) 



1 



Finally, it appears that the best results are obtained when using a n = 9 template. In this case both dE/dt 
and dE/df are over-estimated by the post-Newtonian approximations. However, since /i„ depends on the ratio 
df/dt — {dE / dt) / {dE / df) , it appears that the individual discrepancies are best canceled with 0{v^) expressions. 

The author wishes to apologize for any inconvenience caused by his negligence. This work was supported by the 
Natural Sciences and Engineering Research Council of Canada. 



T. Tanaka, H. Tagoshi, and M. Sasaki, Gravitational wav es by a particle in circular orbit around a Schwarzschild black hole 
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TABLE I. Reduction in signal-to-noise ratio incurred when 
matched filtering with approximate, post-Newtonian tem- 
plates. For each of the considered binary systems, the first 
column lists the order n of the approximation, the second 
column lists TZ„ as calculated using the fully expanded ex- 
pression for Q{v), the third column lists TZ„ as calculated 
using the exact expression for Q(v), and the fourth column 
lists TZn as calculated using the partially expanded expression 
for Q{v). 



n 


Q fully exp 




Q exact 


Q part. exp. 


System A 


(1.4 Mq + 


1.4 Mq) 






4 


0.6756 




0.5852 


0.5217 


5 


0.5565 




0.6117 


0.7157 


6 


0.9578 




0.9451 


0.9320 


7 


0.9777 




0.9866 


0.9802 


8 


0.9731 




0.9772 


0.9782 


9 


0.9848 




0.9854 


0.9858 


10 


0.9836 




0.9837 


0.9838 


11 


0.9836 




0.9837 


0.9838 


System B (1.4 Mq + 


10 Mq): 






4 


0.8580 




0.4890 


0.3772 


5 


0.6718 




0.3788 


0.4305 


6 


0.8271 




0.7455 


0.6859 


7 


0.8938 




0.9977 


0.8796 


8 


0.8762 




0.8945 


0.9111 


9 


0.9962 




0.9650 


0.9459 


10 


0.9897 




0.9958 


0.9984 


11 


0.9924 




0.9977 


0.9995 


System C 


(10 Mq + 


10 Mq): 






4 


0.8367 




0.7610 


0.5855 


5 


0.4433 




0.5354 


0.6684 


6 


0.9769 




0.8500 


0.7988 


7 


0.9038 




0.9997 


0.9469 


8 


0.9006 




0.9314 


0.9541 


9 


0.9998 




0.9850 


0.9669 


10 


0.9941 




0.9983 


0.9995 


11 


0.9959 




0.9992 


0.9999 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9505030v2 
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A particle of mass ^ moves on a circular orbit around a nonrotating black hole of mass M. Under 
the assumption ^ <gi M the gravitational waves emitted by such a binary system can be calculated 
exactly numerically using black-hole perturbation theory. If, further, the particle is slowly moving, 
V = [M^Y^'^ ^ 1 (where v and Q, are respectively the linear and angular velocities in units such that 
G = c — 1) , then the waves can be calculated approximately analytically, and expressed in the form 
of a post-Newtonian expansion. We determine the accuracy of this expansion in a quantitative way 
by calculating the reduction in signal-to-noise ratio incurred when matched filtering the exact signal 
with a nonoptimal, post-Newtonian filter. We find that the reduction is quite severe, approximately 
25%, for systems of a few solar masses, even with a post-Newtonian expansion accurate to fourth 
order — 0{v^) — beyond the quadrupole approximation. Most of this reduction is caused by post- 
Newtonian theory's inability to correctly locate the innermost stable circular orbit, which here is at 
r = 6Af (in Schwarzschild coordinates). Correcting for this yields reductions of only a few percent. 
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Inspiraling compact binary systems, composed of neu- 
tron stars and/or black holes, have been identified 
as the most promising source of gravitational waves for 
interferometric detectors such as LIGO (Laser Interfer- 
ometer Gravitational- wave Observatory [||) and VIRGO 
Q. These systems evolve under radiation reaction, so 
that the gravitational- wave signal increases in amplitude 
as the frequency sweeps through the detector bandwidth, 
from approximately 10 Hz to 1000 Hz. 

Extraction of the information contained in the gravita- 
tional waves, most notably about the masses and spins of 
the companions, will necessitate the construction of accu- 
rate model signals, known as templates 1^-^ . The extrac- 
tion makes use of the well-known technique of matched 
filtering |9|jlC[| , in which the signal is integrated against a 
choice of templates in order to identify the true value of 
the source parameters. It has been established that these 
templates will need to reproduce the signal's phasing es- 
pecially accurately [pPJIlll . This is because the signal 
undergoes a number of oscillations approximately equal 
to 16 000 (if the system is that of two neutron stars) as 
the frequency sweeps through the detector bandwidth. If 
the template loses phase with respect to the true signal, 
even by so little as one cycle, then the signal-to-noise 
ratio is severely reduced, and the information-extraction 
strategy severely impeded. Thus an accuracy of at least 
one part in 16 000 is required for the template's phasing. 

The construction of accurate templates is currently the 
goal of many gravitational- wave theorists |l2) . It is clear 
that the calculation must be based on some approxima- 
tion to the equations of general relativity; the exact inte- 
gration of the equations governing the evolution of com- 
pact binary systems is not currently amenable to numer- 
ical methods. The favored approach is post-Newtonian 
theory, which is based upon an assumption of slow mo- 



tion: If V denotes the orbital velocity, M the total mass, 
and r the orbital separation, then it is assumed that the 
orbital evolution is such that v ^ {M/ry^^ <^ 1. We set 
c = G = 1. 

To date, templates have been calculated accurately 
through second post-Newtonian order p3|-p^, that is, 
order beyond the leading-order, quadrupole-formula 
expression. The question considered in this paper is 
whether this calculation is accurate enough for the pur- 
pose of information extraction. We shall see, as other 
authors have argued |p^-p^, that the answer is nega- 
tive. The post-Newtonian expansion converges extremely 
slowly, if at all. Calculations must be pushed to a much 
higher order. 

The exact treatment of the orbital evolution of a com- 
pact binary system whose companions have comparable 
masses is currently beyond reach. The question posed in 
the preceding paragraph cannot, therefore, be answered 
in this context. However, if instead we consider a system 
composed of a particle with small mass orbiting a black 
hole with much larger mass, then the question can sat- 
isfactorily be answered. This limiting case indeed allows 
for an exact treatment, which can then be compared with 
an approximate, slow-motion treatment. The degree of 
validity of the approximation can thereby be ascertained. 
Our hope is that our conclusions, based upon the small- 
mass-ratio limit, will remain qualitatively valid for sys- 
tems with large mass ratios. 

The gravitational waves generated by a particle in cir- 
cular orbit around a black hole have been the topic of this 
series of papers |l6| , |l9| -p2| . If fj, denotes the particle mass, 
and A'l that of the black hole (which is here assumed to 
be nonrotating), then the assumption fi/M <^ 1 ensures 
that the gravitational perturbations produced by the or- 
biting particle are small and governed by Teukolsky's lin- 
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ear wave equation |g3[ . This equation can be integrated 
using straightforward numerical methods, and the grav- 
itational waveform is thus determined exactly. 

We will focus mainly on the phasing of the waves, 
which is determined by df/dt, the rate at which the 
gravitational-wave frequency changes with time. For cir- 
cular orbits / — r2/7r, where — d(j)/dt is the angular 
velocity, and df / dt can be expressed as 



df _ dE/dt 
dt ^ dE/df 



(1) 



Here, dE/dt is the rate at which gravitational waves re- 
move energy from the system, and dE/df expresses the 
relation between orbital energy and wave frequency. 
The solid curve in Fig. 1 is a plot of 



{dE/dt)N 



where 



1/3 



(2) 



(3) 



is the orbital velocity and {dE/dt)^ = -32ij'^v^° /BW^ 
the Newtonian, quadrupole-formula expression for the 
gravitational-wave luminosity. The plot was obtained 
by numerically integrating the Teukolsky equation, along 
the lines described in Ref. jl6|. An expression for dE/df 
can be obtained by integrating the geodesic equations for 
circular orbits. The orbital energy is defined as _E = ~pt, 
were is the particle's four momentum; it is a constant 
of the motion in the absence of radiation reaction. Wc 
obtain 



Q{v) = 



idE/df)N 



= (1 - 6u2)(l - 3^2) 



-3/2 



(4) 



where {dE/df )n = —tt^M/3v is the Newtonian expres- 
sion. Equation (^) and the solid curve in Fig. 1 give 
through Eq. (Q) an exact representation of df /dt. 

To produce the solid curve in Fig. 1 the Teukolsky 
equation had to be integrated numerically. If, however, 
the small-mass-ratio approximation is combined with a 
slow- motion approximation {v ^ 1), then the Teukol- 
sky equation can be integrated analytically |l^Jl^,|2^,Q . 
The resulting approximate result for P{v) is 

P{v) = 1 - 3.711t;2 -f 12.56v^ - 4.928i;'' - 38.29t>^ 
-I- (115.7 - 16.30 ln^;)w^ - 101. 
-(117.5 - 52. 741nt;)u^ + ---; (5) 

the various coefficients can be found in their analytic 
form in Ref. |jl^. Equation (|^) takes the form of a 
post-Newtonian expansion for P{v). The first five terms, 
through 0{v^), reproduce in the small-mass-ratio limit 
the post-Newtonian calculation of Refs. 25|; the re- 
maining terms have not yet been calculated using post- 
Newtonian theory. A post-Newtonian expansion can also 
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FIG. 1. Various representations of {dE / dt) / {dE / dt) m as 
a function of orbital velocity v = (M/r)^/^ = (ttA//)^''^. 
The solid curve represents the exact result P{v), as cal- 
culated numerically. The various broken curves rep- 
resent the post-Newtonian approximations P^"\v), for 
n — {4, 5, 6, 7, 8}. The smallest value of v corresponds to an 
orbital radius r of 175M; the largest value of v corresponds 
to r = 6M, the innermost stable circular orbit. 



be given for Q{v), by simply expanding Eq. 
V = 0. The result is 



Q{v) = 1 



675, ,6 _ 19845 , ,8 , 



16 



128 



about 



(6) 



Equations (|^) and (^) give through Eq. (^ an approxi- 
mate (post-Newtonian) representation of df /dt. 

Our goal in the sequel is to determine to which ex- 
tent the post-Newtonian representation of df /dt repro- 
duces the true phasing of the waves. We shall denote by 
P{v) and Q{v) the exact version of these functions, while 
Pn(v) and Q„{v) will denote the post-Newtonian repre- 
sentations truncated to the nth power of v. For example, 
Pi{v) represents the right-hand side of Eq. (||) with all 
terms of order and higher removed. 

Figure 1 displays the various curves Pn{v) for n ranging 
from 4 to 8. It is seen that for v > 0.2, these curves 
reproduce rather poorly the exact curve P(v). Moreover, 
the convergence of the post-Newtonian expansion is also 
seen to be poor: adding a term in the expansion does not 
necessarily make the expression more accurate. Witness 
in particular the poor quality of P5{v) compared with 
that of P4{v). [Compare also Ps{v) to Priv).] We shall 
now attempt to determine, in a quantitative manner, how 
much of an obstacle this slow rate of convergence poses 
to information-extraction strategies. 

The starting point of our analysis is an expression for 
the reduction in signal-to-noise ratio incurred when the 
matched filtering of a gravitational- wave signal is carried 
out with a nonoptimal filter. If h{t) denotes the true 
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gravitational waveform, then optimal filtering is achieved 
by using h{t) as a template. This is well-known to yield 
S/N\rna^, the largest possible value of the signal-to-noise 
ratio If, instead, we adopt the post-Newtonian 

approximation hn(t) as a template, then filtering is not 
optimal, and S'/A^j actual, the actual value of the signal-to- 
noise ratio, is smaller than the maximum possible value. 
The reduction in signal-to-noise ratio can be calculated 
to be l26| 



S/N\ 



actual 



mK)\ 



S/N\, 



(7) 



where the inner product (-j-) will be defined presently. 

We pause here to remark that in our calculations, we 
shall not allow the mass parameters /z and M to take dif- 
ferent values in h{t) and hn{t). The problem considered 
here is not that of maximizing the signal-to-noise ratio 
over the source parameters in order to estimate their true 
value. We therefore do not address the issue of the sys- 
tematic errors introduced in parameter estimation when 
using templates which are only approximations to the 
true general-relativistic signal [^. Our considerations 
are less ambitious: the source parameters are assumed 
to be given, and the accuracy of the post-Newtonian 
templates is quantitatively measured by 72.„, the rela- 
tive reduction in signal-to-noise ratio. This measure, we 
believe, is more satisfactory than the mere counting of 
the number of wave cycles contributed by various terms 
in the post-Newtonian expansion of the wave's phasing 

We now return to the definition of the inner prod- 
uct. Let h{f) and hn{f) denote, respectively, the Fourier 
transforms of the functions h{t) and hn{t), with the con- 
vention g{f) ~ J 5(t)e^^*'^* dt for Fourier transforms. Let 
also S{f) be the spectral density of the detector noise, as- 
sumed to be a stationary, Gaussian random process. For 
detectors of the advanced-LIGO type it is appropriate to 
set § 



Sif) = ISo ifo/f) 



Ulkf 



(8) 



for / > 10 Hz; for / < 10 Hz we take S{f) oo. In 
Eq. (^, Sq is a normalization constant irrelevant for our 
purposes, and /o is the frequency at which S{f) is mini- 
mum; we set /o = 70 Hz. The inner product introduced 
in Eq. (Q) is then given by 



{g\h) = 2 



Sif) 



df, 



(9) 



where an asterisk denotes complex conjugation. 

We must now specify h{f) and hn{f), the gravitational 
waveforms in the frequency domain. We use an approx- 
imation [0 in which the amplitude of the waveform is 
described accurately to Newtonian order, while its phase 
is described exactly in the case of h{f), and by a post- 
Newtonian expansion in the case of hn{f). Consequently, 
we have gf] 



/i(/)-^r'/'exp[jV'(/)], 
^n(/)=^/"'/'exp[zV'„(/)], 



(10) 



where ^ is a constant. The phase functions il){f) and 
i'nif) are constructed as follows. 

We consider ipif)] '4'n{f) is dealt with similarly. Our 
starting point is the wave's phasing in the time domain, 
which is determined by df /dt given above. Combining 
Eqs. (|l|)-(|) we obtain df/dt = Q&^iv^^ P{v) /hirM^Qiv). 
Integration then yields 



t{v)lM = tjM + — 



QW) 



dv' 



for time as a function of velocity, and 
5M r QW) 



dv' 



(11) 



(12) 



for the phase = / 2'Kfdt. In Eqs. (O) and (p^), Vi is an 
arbitrary reference point, andt^ = t(ui), $i = $(ui). The 
frequency-domain phase function ^{f) is obtained via 
the stationary phase approximation [^J^,|28t, according 



to which ?/)(/) = 2-Kft{v) - $(«) - 7r/4. Wc therefore find 
V'(/) = 2(t,/M)«3 _ _ 7r/4 



5M r [v^^ - v"^)Q{v') 
16^ X. t)'9P(w') 



dv'. 



(13) 



where, we recall, v = (ttM/)^/"^. Equations ( pi] ) and 
(p^), together with the fact that the gravitational- wave 
signal must be cut off at a frequency /isco correspond- 
ing to the innermost stable circular orbit, completely 
specify the Fourier transform of the waveform. We use 



(M/r-isco)^/^ = 6-3/2, ^.^ 



6Af is 



the radius of the innermost stable circular orbit of the 
Schwarzschild spacetime. 

The calculation of TZn can now be carried out. 
Straightforward manipulations, using Eqs. (|7|)-(p^), yield 



1/7 



X '^/'^ COS At/) 
+ 2 + 2x2 



where x = f/fa, xisco = /isco//o, and 

A?/' = V' - V'n- 

The constant / in Eq. ( p4[ ) is given by 

2.-7/3 



1/7 



2x2 



■ dx. 



(14) 



(15) 



(16) 



and ensures that TZn = 1 if '0(/) = V'n(/)- 

The numerical value of 7^„, for given n, depends on 
the value of the constants ti and $i which appear in 
4'{f)i and on the value of the constants i„i and $„i which 
appear in ipnif)- To maximize TZn, we set ti ~ tni and 
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choose — such that Atfj vanishes at the value of 
/ for which x~'^^^{x~'^ + 2 + 2x^)~^ is maximum. It is 
easy to check that this occurs at x^^ax — 0.6654, so that 
/max — 46.58 Hz. A simple calculation also shows that 
W = (7rM/max)i/3 =^ O.O8966(M/Af0)i/3, where M© 
denotes the mass of the Sun. With these choices, Aip 
becomes 



^ , 5M r v^-v'' 



Q{V') Qnjv') 
P{V') Pn{v') 



dv'. (17) 



It is a straightforward numerical problem |29| to compute 
At/j for a given v = {-kM fY^^ and to then evaluate the 
integral to the right-hand side of Eq. (p^. 

The calculation presented above is, strictly speaking, 
only applicable to binary systems with small mass ra- 
tios. In this limit, the only place in which fi/M appears 
is as an overall multiplicative factor in Eq. (pT[). In this 
limit, therefore, the phase lag Aij} scales exactly as M/ ji. 
As the mass ratio is allowed to increase, our results for 
P{v), Q{v), and their post-Newtonian analogues, must 
be corrected p8L For example, the constant coefficients 
in Eq s. (Jq ) and(|) would acquire /x/M-dependent correc- 
tions | |13[ |. Our expressions for {dE/dt)N and {dE/df)N, 
however, stay valid for large mass ratios, provided that 
/i is then interpreted as the system's reduced mass, and 
M as the system's total mass. 

In the following we will let n/M become large, with- 
out modifying our expressions for P{v), Q{v), and their 
post-Newtonian analogues. We will assume, without jus- 
tification, that the qualitative behavior of these functions 
is not appreciably affected by the finite-mass-ratio cor- 
rections [which are not known, apart from the terms of 
lowest order in P„(w) and Qn{v)]- We will therefore apply 
our formalism to binary systems with comparable masses, 
in the hope that our conclusions based on the /i/M ^ 
limit will be qualitatively valid also in the large-mass- 
ratio case. We note, in accordance with our previous 
observation, that in Eq. ( |l7|) M is to be interpreted as 
the total mass and /i as the reduced mass. 

Our results are summarized in Table I. We consider 
three types of binary systems. The first (system A) con- 
sists of two neutron stars, each with a mass equal to 
1.4M0. The second (system B) consists of one neutron 
star (I.4M0) and one black hole (IOMq). The third (sys- 
tem C) consists of two black holes (lOM© each). 

The second column in Table I lists TZn as defined in 
Eqs. (0) and (0, for n ranging from 4 to 8. This shows 
that for system A, the signal-to-noise ratio is always 
smaller than 0.7651 times the maximum possible value, 
with the best result obtained when n = 7, which cor- 
responds to dE/dt accurate to 3.5PN order, and dE/df 
accurate to 3PN order. The reduction in signal-to-noise 
ratio is therefore quite severe, even at such a high order 
in the post-Newtonian expansion. The corresponding re- 
sults for systems B and C can also be obtained from Table 
I. 



TABLE I. Reduction in signal-to-noise ratio incurred when 
matched filtering with approximate, post-Newtonian tem- 
plates. For each of the considered binary systems, the first 
column lists the order n of the approximation, the second 
column lists T?.'"' as calculated using the post-Newtonian ap- 
proximation Q'"-* for dE/df, the third column lists T?.*-"' as 
calculated using the exact expression Q{v) for dE/df, and the 
fourth column lists T?.'"' as calculated using the alternative 
post-Newtonian approximation Q'^"'\v) for dE/df. 



n '(v) 






Q (v) 


bystem A Mq 


+ 


i.4 Mq): 




4 0.5796 




0.4958 


0.4326 


5 0.4646 




0.5286 


0.6492 


6 0.7553 




0.9454 


0.9300 


7 0.7651 




0.9864 


0.9819 


8 0.7568 




0.9695 


0.9709 


System B (1.4 M© 


+ 


10 Mq): 




4 0.8478 




0.3954 


0.2916 


5 0.2413 




0.2922 


0.4012 


6 0.6097 




0.6788 


0.6051 


7 0.6023 




0.9966 


0.8528 


8 0.5734 




0.8744 


0.8942 


System C (10 Mq 


+ 


10 Mq): 




4 0.8232 




0.6919 


0.4910 


5 0.3515 




0.4422 


0.5866 


6 0.7107 




0.8088 


0.7417 


7 0.7201 




0.9997 


0.9310 


8 0.6730 




0.9152 


0.9420 



The reduction in signal-to-noise ratio is due to the fact 
that Pn{v) is only an approximation to P{v), see Fig. 1, 
and that Qn{v) is only an approximation to Q(v). It is 
interesting to ask how much of the signal-to-noise ratio 
could be recovered if only P{v) were approximated by a 
post-Newtonian expansion, while Q{v) were kept exact. 
To answer this question amounts to repeating the cal- 
culation presented previo usly , but with Q{v) substituted 
in place of Qn{v) in Eq. (p^. The results are shown in 
the third column of Table I. Not surprisingly, we see that 
keeping dE/df exact gives much better results. For sys- 
tem A, TZn can become as large as 0.9864 (for n — 7), and 
is already as large as 0.9454 for n — 6 (which corresponds 
to dE/dt accurate to 3PN order). The corresponding re- 
sults for systems B and C can also be obtained from Table 
I. 

Why does the exact expression for dE/df give such 
better results? The answer is that while the exact expres- 
sion for dE/df correctly vanishes at v = Uisco = 6^^^^ (at 
the innermost stable circular orbit), its post-Newtonian 
analogue fails to do so. For example, Qs{v) goes to 
zero at w ~ 0.4236 > 6^^/^, corresponding to a radius 
r ~ 5.572M. 

To establish that this indeed the reason, we consider 
an alternative expression for Q{v), obtained from Eq. (^) 
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by expanding only the (1 — 3u^) factor: 

Q'(z;) = (l-6«2)(l + |«2+ 135^4 



945 „. 6 , 25515 „,8 
16 '■^ 128 ^ 



(18) 



This expression, together with its truncated versions 
Q'j^{v), manifestly go to zero at u = 6"^/^. We have re- 
peated the calculations with Q'n{v) substituted in place 
of Qni^) in Eq. (17). The results are shown in the fourth 
column of Table I. It is evident that, generally speaking, 
using Q'n(v) gives much larger values of TZn than using 
Qn{v). This is especially remarkable for system A, since 
the relative contribution to the signal-to-noise ratio com- 
ing from large values of v is very small: For v = 
and M = 2.8Af0 the corresponding frequency is 1570 Hz, 
at which S{f ) is 200 times its minimum value Sq. 

These results suggest that accurate knowledge of the 
location of the innermost stable circular orbit (where 
dE/df goes to zero) might significantly improve the per- 
formance of the post-Newtonian templates, through a 
factorization of the kind shown in Eq. ( [T8| ) . It is possible 
that in the case of binary systems with large mass ratios, 
such information could be obtained by numerically solv- 
ing the initial value problem of general relativity |pO| , ^ , 
which would be far less laborious than solving the full 
dynamical problem. Another approach would be to use 
the Kidder- Will- Wiseman "hybrid" equations of motion 
p2| , which give a better approximation to the innermost 
stable circular orbit than the standard post-Newtonian 
equations. 

Conversations with Cliff Will were greatly appreci- 
ated. This work was supported by the Natural Sci- 
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